A Process Calculus for Molecular Interaction Maps

We present the MIM calculus, a modeling formalism with a strong biological basis, which provides biologically-meaningful operators for representing the interaction capabilities of molecular species. The operators of the calculus are inspired by the reaction symbols used in Molecular Interaction Maps (MIMs), a diagrammatic notation used by biologists. Models of the calculus can be easily derived from MIM diagrams, for which an unambiguous and executable interpretation is thus obtained. We give a formal definition of the syntax and semantics of the MIM calculus, and we study properties of the formalism. A case study is also presented to show the use of the calculus for modeling biomolecular networks.Comment: 15 pages; 8 figures; To be published on EPTCS, proceedings of MeCBIC 200

Depicting combinatorial complexity with the molecular interaction map notation

To help us understand how bioregulatory networks operate, we need a standard notation for diagrams analogous to electronic circuit diagrams. Such diagrams must surmount the difficulties posed by complex patterns of protein modifications and multiprotein complexes. To meet that challenge, we have designed the molecular interaction map (MIM) notation (http://discover.nci.nih.gov/mim/). Here we show the advantages of the MIM notation for three important types of diagrams: (1) explicit diagrams that define specific pathway models for computer simulation; (2) heuristic maps that organize the available information about molecular interactions and encompass the possible processes or pathways; and (3) diagrams of combinatorially complex models. We focus on signaling from the epidermal growth factor receptor family (EGFR, ErbB), a network that reflects the major challenges of representing in a compact manner the combinatorial complexity of multimolecular complexes. By comparing MIMs with other diagrams of this network that have recently been published, we show the utility of the MIM notation. These comparisons may help cell and systems biologists adopt a graphical language that is unambiguous and generally understood

A formal MIM specification and tools for the common exchange of MIM diagrams: an XML-Based format, an API, and a validation method

<p>Abstract</p> <p>Background</p> <p>The Molecular Interaction Map (MIM) notation offers a standard set of symbols and rules on their usage for the depiction of cellular signaling network diagrams. Such diagrams are essential for disseminating biological information in a concise manner. A lack of software tools for the notation restricts wider usage of the notation. Development of software is facilitated by a more detailed specification regarding software requirements than has previously existed for the MIM notation.</p> <p>Results</p> <p>A formal implementation of the MIM notation was developed based on a core set of previously defined glyphs. This implementation provides a detailed specification of the properties of the elements of the MIM notation. Building upon this specification, a machine-readable format is provided as a standardized mechanism for the storage and exchange of MIM diagrams. This new format is accompanied by a Java-based application programming interface to help software developers to integrate MIM support into software projects. A validation mechanism is also provided to determine whether MIM datasets are in accordance with syntax rules provided by the new specification.</p> <p>Conclusions</p> <p>The work presented here provides key foundational components to promote software development for the MIM notation. These components will speed up the development of interoperable tools supporting the MIM notation and will aid in the translation of data stored in MIM diagrams to other standardized formats. Several projects utilizing this implementation of the notation are outlined herein. The MIM specification is available as an additional file to this publication. Source code, libraries, documentation, and examples are available at <url>http://discover.nci.nih.gov/mim</url>.</p

The NCI-60 methylome and its integration into CellMiner

A unique resource for systems pharmacology and genomic studies is the NCI-60 cancer cell line panel, which provides data for the largest publicly available library of compounds with cytotoxic activity (∼21,000 compounds), including 108 FDA-approved and 70 clinical trial drugs as well as genomic data, including whole-exome sequencing, gene and miRNA transcripts, DNA copy number, and protein levels. Here, we provide the first readily usable genome-wide DNA methylation database for the NCI-60, including 485,577 probes from the Infinium HumanMethylation450k BeadChip array, which yielded DNA methylation signatures for 17,559 genes integrated into our open access CellMiner version 2.0 (https://discover.nci.nih.gov/cellminer). Among new insights, transcript versus DNA methylation correlations revealed the epithelial/mesenchymal gene functional category as being influenced most heavily by methylation. DNA methylation and copy number integration with transcript levels yielded an assessment of their relative influence for 15,798 genes, including tumor suppressor, mitochondrial, and mismatch repair genes. Four forms of molecular data were combined, providing rationale for microsatellite instability for 8 of the 9 cell lines in which it occurred. Individual cell line analyses showed global methylome patterns with overall methylation levels ranging from 17% to 84%. A six-gene model, including PARP1, EP300, KDM5C, SMARCB1, and UHRF1 matched this pattern. In addition, promoter methylation of two translationally relevant genes, Schlafen 11 (SLFN11) and methylguanine methyltransferase (MGMT), served as indicators of therapeutic resistance or susceptibility, respectively. Overall, our database provides a resource of pharmacologic data that can reinforce known therapeutic strategies and identify novel drugs and drug targets across multiple cancer type

Molecular Dynamics Simulations

A tutorial introduction to the technique of Molecular Dynamics (MD) is given, and some characteristic examples of applications are described. The purpose and scope of these simulations and the relation to other simulation methods is discussed, and the basic MD algorithms are described. The sampling of intensive variables (temperature T, pressure p) in runs carried out in the microcanonical (NVE) ensemble (N= particle number, V = volume, E = energy) is discussed, as well as the realization of other ensembles (e.g. the NVT ensemble). For a typical application example, molten SiO2, the estimation of various transport coefficients (self-diffusion constants, viscosity, thermal conductivity) is discussed. As an example of Non-Equilibrium Molecular Dynamics (NEMD), a study of a glass-forming polymer melt under shear is mentioned.Comment: 38 pages, 11 figures, to appear in J. Phys.: Condens. Matte

Density functional method for nonequilibrium electron transport

We describe an ab initio method for calculating the electronic structure, electronic transport, and forces acting on the atoms, for atomic scale systems connected to semi-infinite electrodes and with an applied voltage bias. Our method is based on the density functional theory (DFT) as implemented in the well tested Siesta approach (which uses non-local norm-conserving pseudopotentials to describe the effect of the core electrons, and linear combination of finite-range numerical atomic orbitals to describe the valence states). We fully deal with the atomistic structure of the whole system, treating both the contact and the electrodes on the same footing. The effect of the finite bias (including selfconsistency and the solution of the electrostatic problem) is taken into account using nonequilibrium Green's functions. We relate the nonequilibrium Green's function expressions to the more transparent scheme involving the scattering states. As an illustration, the method is applied to three systems where we are able to compare our results to earlier ab initio DFT calculations or experiments, and we point out differences between this method and existing schemes. The systems considered are: (1) single atom carbon wires connected to aluminum electrodes with extended or finite cross section, (2) single atom gold wires, and finally (3) large carbon nanotube systems with point defects.Comment: 18 pages, 23 figure

Predicted Functions of MdmX in Fine-Tuning the Response of p53 to DNA Damage

Tumor suppressor protein p53 is regulated by two structurally homologous proteins, Mdm2 and MdmX. In contrast to Mdm2, MdmX lacks ubiquitin ligase activity. Although the essential interactions of MdmX are known, it is not clear how they function to regulate p53. The regulation of tumor suppressor p53 by Mdm2 and MdmX in response to DNA damage was investigated by mathematical modeling of a simplified network. The simplified network model was derived from a detailed molecular interaction map (MIM) that exhibited four coherent DNA damage response pathways. The results suggest that MdmX may amplify or stabilize DNA damage-induced p53 responses via non-enzymatic interactions. Transient effects of MdmX are mediated by reservoirs of p53∶MdmX and Mdm2∶MdmX heterodimers, with MdmX buffering the concentrations of p53 and/or Mdm2. A survey of kinetic parameter space disclosed regions of switch-like behavior stemming from such reservoir-based transients. During an early response to DNA damage, MdmX positively or negatively regulated p53 activity, depending on the level of Mdm2; this led to amplification of p53 activity and switch-like response. During a late response to DNA damage, MdmX could dampen oscillations of p53 activity. A possible role of MdmX may be to dampen such oscillations that otherwise could produce erratic cell behavior. Our study suggests how MdmX may participate in the response of p53 to DNA damage either by increasing dependency of p53 on Mdm2 or by dampening oscillations of p53 activity and presents a model for experimental investigation

Functional Categories Associated with Clusters of Genes That Are Co-Expressed across the NCI-60 Cancer Cell Lines

The NCI-60 is a panel of 60 diverse human cancer cell lines used by the U.S. National Cancer Institute to screen compounds for anticancer activity. In the current study, gene expression levels from five platforms were integrated to yield a single composite transcriptome profile. The comprehensive and reliable nature of that dataset allows us to study gene co-expression across cancer cell lines.Hierarchical clustering revealed numerous clusters of genes in which the genes co-vary across the NCI-60. To determine functional categorization associated with each cluster, we used the Gene Ontology (GO) Consortium database and the GoMiner tool. GO maps genes to hierarchically-organized biological process categories. GoMiner can leverage GO to perform ontological analyses of gene expression studies, generating a list of significant functional categories.GoMiner analysis revealed many clusters of coregulated genes that are associated with functional groupings of GO biological process categories. Notably, those categories arising from coherent co-expression groupings reflect cancer-related themes such as adhesion, cell migration, RNA splicing, immune response and signal transduction. Thus, these clusters demonstrate transcriptional coregulation of functionally-related genes